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The problem of the one-exciton absorption spectrum is considered for the lattice 
of two-level interacting atoms whose initial energy splitting depends on the coordi- 
nate. It is shown that for some types of interatomic interaction, this problem can be 
reduced to a differential equation of Schrodinger type, which, in some cases, can be 
solved in closed form. By way of example, problems on the one dimensional chain 
having an initial splitting jump and on the spherical cluster are solved. In both cases 
analitical solutions are obtained, which considerably reduce the time required for cal- 
culations of one-dimensional systems and permit the calculation of three-dimensional 
spherical clusters, whereas the numerical calculation becomes impossible for the clus- 
ters with a radius larger than ten lattice constants. The spectra calculated for small 
systems completely coincide with the spectra obtained by means of numerical diago- 
nalization. The proposed method may appear to be convinient for solving problems 
related to Frenkel excitons in systems having no translation symmetry (disordered 
systems) . 

1. Introduction 

The recent advent of technology of manufacturing small-sized crystalline structures has stim- 
ulated spectroscopic studies of such objects. These objects include, for example, quantum 
wells, quantum wires, J-aggregates and nanocrystals. Physical properties of such objects 
(absorption spectra, the specific heat, etc.) are significantly determined by their restricted 
sizes. From the point of viev of optical spectroscopy, such systems can be treated, in some 
cases, as fragments of crystal lattice, and their absorption spectra can be interpreted using 
the concept of excitons. In this case it is useful to have simple computational technique 
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for estimating the scale of effects related to the restricted size of these systems. This paper 
deals with one such method. 

In the second section, the mathematical scheme of the method is presented, and the 
differential equation is obtained for the exciton Green's function required for calculation 
of the absorption spectrum. In the third section, the one dimensional exponential model 
of an exciton is formulated, in which the interatomic interaction exponentially varies in 
space, and the equation obtained in the second section has a simpler form and can be 
exactly solved in some cases. As an example, expressions are obtained for the absorption 
spectra of a finite chain and a chain in which the distribution of the initial splitting energy 
has the form of meander. In the forth section, the suggested method is generahzed to 
the case of three-dimensional restricted systems. In this case, the equation for Green's 
function simplified, if the interatomic interaction is discribed by the Yukawa potential. In 
this section, the expression derived for the absorption spectrum of a spherical cluster. The 
spectra obtained in Sections 3 and 4 exhibit a well developed structure and strongly differ 
from the monochromatic spectrum of an infinite crystal. Almost all spectra were tested 
by means of a direct computer diagonalization, and, in all cases, complete agreement was 
achieved. The computer calculation of the spectra of one-dimensional chains with lenths 
exceeding 100 lattice constants a takes much more time than the calculation according the 
above expressions. Numerical computer-aided analysis of the spectral structure of a three- 
dimensional cluster requires several hours for the clusters radius 10a. At the same time, as 
shown in Section 4, the spectrum of a cluster with radius of 30a exhibits the well-developed 
structure related to the restricted size of the cluster. The expressions obtained in the Section 
4 can be used for estimating effects caused by the restricted size of aspherical clusters. 

Note that all the expressions obtained in this paper are essentially based on the expo- 
nential (one-dimensional problems) and Yukawa (three-dimensional problem) types of inter- 
atomic interaction. This, however, does not appreciably restrict the field of application of the 
above method for the following reason. The exciton spectrum is often discribed in effective 
mass approximation, which adequately reflects the dispersion law only for long-wavelength 
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excitons. The presence of two parameters in exponential (Yukawa) potential allows one to 
simulate not only the effective mass of an exciton but also the exciton bandgap. 



2. Equation for the exciton Green's function 

Consider a regular lattice of two-level atoms. We assume that the initial (neglecting the 
interatomic interaction) splitting of atomic levels is a given function e{r) of the atomic co- 
ordinate r. Consider also the interatomic interaction, which can cause the energy transfer 
between atoms, the dependence of this interaction on the interatomic distance being deter- 
mined by the function w{r) We will study only one particle excited states of this system 
and the absorption spectrum ralated to the transitions to these states from the ground 
state, i.e. the one-exciton absorption spectrum. If the transition dipole moment is the same 
for all atoms in the lattice, this spectrum can be found by diagonalization of the matrix 
Hrr' = Srr'£{r) + w{r — r'), the positions of spectral lines being determined by eigenvalues, 
and their intensites-by the squares of sums of components of corresponding column eignvec- 
tors (it is in this way that squares of moduli of corresponding matrix elements upon transi- 
tions from ground state are calculated). Below, however, we will use the following method 
of calculation of the above spectrum. We will construct Green's function grr' = {E — H)~^, 
for matrix H and define the function G{r, E) as follows 

G{r,E) = J2gri 

Then, the absorption spectrum A{E) , correct to insignificant factors , has the formB 
A{E) = -Im J2 G{r, E + i6) (1) 

r 

each line in the spectrum being discribed by Lorentzian with the half-width 6. By 
considering w{r — r') in the matrix if as a perturbation, we write Dyson equation for 
Green's function 
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9rr' 



I 

Summation over r' yields the equation for the function G in which we will pass from 
summation over I to the integration: 



Here p is the lattice density. The passage from summation to integration in (2) is the 
only approximation we use, and it should be noted that such passage is acceptable. The 
sum over I in the Dyson equation is close to the integral in (2), provided the function 
w{r — r')G{r') changes only slightly when r' changes by a lattice constant a. For this to 
happen, first, it is required that R > a. Second, G{r) also should not strongly change 
on a scale ~ a, which should be checked after calculations according to the procedure 
suggested below. Note that in the region of absorption spectrum, the function G oscillates 
at the spacial frequency of the resonant exciton. Because the main contribution to the 
absorption spectrum (1) is made by long- wavelength excitons, the passage from summation 
to integration in (2) is justified. It should be emphasized that the function G in (1) describes 
the absorption spectrum rather than the density of states whose calculation would require 
the consideration of short- wavelenght excitons. 

Let us now define the function ^(r) = G{r){E — e{r)) It follows from (2) that 



Integral equation (3) can be reduced to the differential equation in the following way. 
One can construct the differential operator L{d/dr) for wide class of functions w{r), so that 
w{r) will be the source function (Green's function) for this operator, i.e.. 




(2) 




(3) 



where 
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L{d/dr)w{r) = S{r) 

To do this, equation (4) should be written in the Fourier representation 



(4) 



w{r) = — I e^^''w{k)dk, 



L{jk)w{k) = 1 



(5) 



The function L is determined from the latter equation. Consider now the inhomogeneous 
equation 



L{d/dr)^ = $(r) 



(6) 



The particular solution of this equation can be obtained with the help of the source 
function w(r): 



^'(r) = j w{r - r')^{r')dr' 



(7) 



If 



E - e{r\ 

then expression (7) coincides with (3) . Therefore, the required solution ^{r) of integral 
equation (3) satisfies the differential equation 

\l/(r) 



L(a/9r)*(r) = p 



E-eir) 



(8) 



In a number of cases, it is easier to solve the equation (8) than (3). If the function is 
known, the absorption spectrum is detrmined by an expression that is similar to (1): 



A{E) = -Im j 



E + jd-e{r 



-dr 



(9) 



where S is lattice area. 
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3. One-dimensional chain with exponential interatomic interaction 



Consider a one-dimensional chain with interatomic interaction of the form w{r) — 
V exp{—\r\/R). One can obtain from (5) the following expression for the operator L{d/dr) 



L{d/dr) 



E_ 

W 



R 



-2 



d''/dr' 



and equation (8) is transformed to the equation 



(10) 



(5^^R 



W 



- 1 



E - e{r) 
W = 2pVR, 

Which is similar to the inhomogeneous Schrodinger equation with the 'potential' /3^(r). 
For this reason, the function ^ has a continuous first derivative. Let us find \1' for a finite 
chain whose ends have coordinates ±1. In this case, e{r) has the form 



e{r) 



oo r < —I 
oo r > I 



(.0 -I <r <l 

i.e., the region outside the chain is as if filled by infinitely detuned atoms to which the 
excitation cannot arrive. Of course, this is only the 'variable' component of the splitting. 
The 'constant' component results in a shift of the spectrum as a whole, and we will not take 
it into account below. In each region , the solution of (10) can be obtained quite easily. 
Taking into account the symmetry of the chain, the solution can be written in the form 



' Ci cos(/?r) - l/(/3i?)2 -l<r <l 
C2exp(-r/i?) + 1 r>l 
,C2exp(r/i?) + l r<-l 



where 
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If E is not too close to zero (see below) , then /3 is the wave number of an exciton with 
energy E. In regions outside the chain, we retained only decreasing exponentials, because 
the function G vanishes there. This can be done passing to the limit £ — > oo for r < — / and 
r > Z, if the function ^ is limited. Constants Ci and C2 are determined by the continuity 
condition for the function ^ and its derivative at the chains ends. By dropping simple 
calculations, we present the final result for the function G inside the chain: 

air) = ^(..^ ^,^?:ri ^ (U) 



E-W\ E pi sm{pl) - cos{pl) . 
Let us indicate the region of energies E where the absorption spectrum calculated from 

(1) and (11) can be used. The passage from summation to integration in Dyson equation 

(2) will be justified, if the spacial half-period of G (equal to vr//?) is grater than the lattice 
constant a. This results in the condition 

Figure 1 (There are no figures in this version, see figures in Optics and Spectroscopy, 
Vol.82, No. 2, 1997, pp 242-246) shows the absorption spectrum calculated from (11) and 
(1) for the chain parameters / = 50, R = 2,V = 1, and 6 = 0.1, the energy E being plotted in 
units of V. The vertical bar along the abscissa axis shows the energy Ec- One can see that 
in the "forbidden" region, the spectrum is in fact absent, and this region is of no interest. 
The spectrum of this chain was also obtained by numerical diagonalization and proved to 
be completely coincident with the spectrum in Fig.l. 

Let us also consider here the somewhat less trivial case of a chain in which e{r) has the 
form of a meander: 

'00 r < —h or r > h 
e{r) = \ e I < r < h or -I > r > -h 
.0 -I <r <l 



We assume for the simplicity of calculations, that this chain is closed to a ring, i.e., the 
point h coincides with —h. The solution of (10) is obtained in similar manner, the only 
difference being that it is necessary to require the values of the function "if and its derivative 
to be the same in the points h and —h. The result for integrated Green's function, whose 
imaginary part determines the absorption spectrum, is 

^fcdr- ^ I f3f-f3l{W-E)-'-{W + e-E)-' 

2 J E-W E-W-e (3i(32 /32ctg(/3i/) + /3ictg(/32(/i - 0) 

—h 

where 

£l,2 = 0,£ 

In this case, condition (12) should be supplemented by the requirement 
e-Ec>E>e + Ec (14) 

because for energies lying within the Ec vicinity of e, the function G has a spatial period 
smaller than a on the chain fragments with initial splitting e. Note, however, that in this 
case, the contribution from the above fragments will be small compared to that from the rest 
of the chain, where the function G oscillates not so rapidly. For this reason, the absorption 
spectrum can be realistic even when condition (4) is violated. 

The spectrum of the chain with a jump of the initial splitting is presented in Fig. 2. The 
parameters arc: h = 30,1 = 10, R = 2,V = 1,6 = 0.1, and e = 2.5; the energy is plotted in 
units of V. The spectrum in Fig. 2a is obtained by numerical diagonalisation, and that in 
Fig 2b is calculated from (13). The regions of forbidden energies are indecated by a vertical 
bar Eq and horizontal bar [condition (14)]. For the chains of lenght ~ 100a, calculation by 
means of the above equations takes much less time than numerical diagonalization. In the 
next section the problem of three-demensional spherical cluster will be solved. The numerical 
calculation of such a claster with radius > 20a represents a vast computer problem, while 
clusters with radii ~ 150a still exhibit the distinct structure, related to their finite size. 
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4. Spherical cluster with interatomic interaction of Yukawa type 

The approach presented in section 2 can be generahzed to the case of an arbitrary de- 
mensionahty. In this section we consider a three-dimensional lattice with an inteaction of 
w{r) = V exp{—r/R)/r (the Yukawa potential). It is well known [and can be easily verified 
from (5)] that 

After calculating the integral A^^ in (3) in the spherical coordinate system, we find that 
in this case, equation (8) again has the form of the inhomogeneous Shredinger equation 

where 





|r| 


[ oo 


r 



and, hence, function ^ is continuous with its first derivative. 

Let us now solve the problem on a spherical cluster with radius cr, i.e., the function e{r) 
will be taken in the form 

r| < (J 

£{r) 

By representing the Laplace operator in (15) in spherical coordinates, for function (it 
has only the radial dependence), we obtain the equation 

where 
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In our case, 



\ {W/E - r<a 
The general solution of equation (16) for both regions has the form 



/ 1 \ ^ sin Br , cos Br 

* = -fe) +"^+'^ 

Taking into account the boundedness of ^ in the region r > cr, we obtain that in this 
region, 

* = l + ai— ^-^ — '—^ (17) 
r 

Function should be also bounded at r = 0. Taking this into account, we obtain that in 
the region r < cr, 

^^EG^ + a,'^ (18) 
r 

where 



G^ = {E-W)-' 

Function EGoc is the solution of (15) in the case of an infinite homogeneous lattice. For 
this reason, lmGoo{E + jS) determines the absorption spectrum of such a lattice. Constants 
Oi and a2 are determined from continuity conditions for the function ^ and its first derivative 
for r = a The absorption spectrum is determined from (9) 



^ ^ J E + j6 



which yields the result 



A(E) = -Imjf^ - W{R + a)[sin aa-aa cos aa] ^a^^^ + jS) \ (19) 
^ ^ IV3 {E + J 5) {Ra cos aa + sin aa)' ^ ' ^ ^ 

where 
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1 . w 

E + j5 



- 1 



Expression (19) contains two contributions: the bulk contribution (the first term in 
parentheses) and the surface contribution (the second term) . The first term gives the spectral 
line that corresponds to an infinite lattice, while the second one is related to the structure 
caused by the finite size of the cluster. 

Figure 3 shows the absorption spectrum of a spherical cluster calculated from (19) for 
the parameters V = 1, R = l,a = 5, and S = 0.15; the energy is plotted in units of V. 
The vertical bar shows the energy Eq- The spectrum obtained by numerical diagonalization 
coincides completely with the spectrum in Fig. 3. However, it should be noted that the 
Yukawa potetial has a singularity at r = 0. For this reason, the question arises upon 
numerical calculations as to what form the diagonal elements of the H^r' matrix should take 
(see Section 2). Although these matrix elements do not affect the shape of the absorption 
spectrum, they may cause its shift. In the numerical verification of the spectrum in Fig. 3, 
we assumed that Hrr' — 0, because in the calculation of intergrals with Yukawa potential, 
the contribution from the region r ~ is insignificant due to the smallness of the volume 
r^dr, and, in fact, the "effective potential" w{r)r^ is used evrywhere, which vanishes for 
r — 0. Computer calculations performed with other values of parameters are also in complete 
agreement with (19). 

Figure 4 shows the spectrum of spherical cluster for V = 1, R = 2,a = 30, and 5 = 
0.6. Comparison with numerical diagonalization was not performed, because, in this case, 
computer calculations are time consuming. Figure 4 shows that despite large radius of the 
cluster, its spectrum differs from that of infinite lattice (the singlet at E = W) and exhibits 
a structure. Clusters with a radius of 150a also show a noticeable structure (on the scale of 
the bandgap W). 
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